Cooperative dynamics of microtubule ensembles: Polymerization forces and 

rescue-induced oscillations 
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We investigate the cooperative dynamics of an ensemble of N microtubules growing against an 
elastic barrier. Microtubules undergo so-called catastrophes, which are abrupt stochastic transitions 
from a growing to a shrinking state, and rescues, which are transitions back to the growing state. 
Microtubules can exert pushing or polymerization forces on an obstacle, such as an elastic barrier if 
the growing end is in contact with the obstacle. We use dynamical mean-field theory and stochastic 
simulations to analyze a model where each microtubule undergoes catastrophes and rescues and 
where microtubules interact by force sharing. For zero rescue rate, cooperative growth terminates 
in a collective catastrophe. The maximal polymerization force before catastrophes grows linearly 
with N for small N or a stiff elastic barrier, in agreement with available experimental results, whereas 
it crosses over to a logarithmic dependence for larger iV or a soft elastic barrier. For a nonzero rescue 
rate and a soft elastic barrier, the dynamics becomes oscillatory with both collective catastrophe and 
rescue events, which are part of a robust limit cycle. Both the average and maximal polymerization 
forces then grow linearly with N, and we investigate their dependence on tubulin on-rates and 
rescue rates, which can be involved in cellular regulation mechanisms. We further investigate the 
robustness of the collective catastrophe and rescue oscillations with respect to different catastrophe 
models. 

PACS numbers: 87.16.Ka, 87.16.-b 



I. INTRODUCTION 

Microtubules (MTs) are long and stiff filamentous 
proteins, which assemble and disassemble from tubulin 
dimers and serve various functions in the cytoskclcton: 
MT stiffness plays an important role in cytoskeletal me- 
chanics but they also serve as "tracks" for intracellular 
transport by molecular motors. Finally, MTs exhibit an 
unusual polymerization dynamics, which is essential for 
spatial organization and remodelling processes in the cy- 
toskeleton 

Polymerizing MTs can also generate pushing forces 
within the cell, because MTs grow by tubulin insertion 
also in the presence of an obstacle which exerts an op- 
posing force on tubulin dimers inserted at the growing 
MT tip. The opposing force is transmitted onto the MT, 
as has been demonstrated in singlc-MT buckling experi- 
ments in front of a solid wall Q . The force slows down 
further MT growth because monomer insertion against 
a force involves additional mechanical work. Growth fi- 
nally terminates under a maximal polymerization force, 
which is typically in the pN range [2|. 

An important and unique feature of the MT poly- 
merization dynamics is their so-called dynamic instabil- 
ity. The dynamic instability gives rise to phases of fast 
shrinking, which stochastically interrupt polymerization 
phases |3( ■ Each phase of of fast shrinking is initiated by 
a catastrophe event and terminated by a rescue event. 
This complex dynamic behavior is central to rapid re- 
modelling in the cytoskclcton but also affects the poly- 
merization force. 

Catastrophe and rescue events of MTs are associated 
with guanosinc triphosphate (GTP) hydrolysis within the 



MT Q. Each tubulin dimer consists of an a- and /3- 
tubulin and contains two GTP binding sites. Tubulin 
dimers assemble into polar MTs with a fast-growing plus 
end. GTP-tubulin dimers attach to the plus end with 
both binding sites containing GTP. The GTP in the a- 
tubulin is hydrolyzed to guanosine diphosphate (GDP) 
within the MT. This process leads to the formation of a 
GTP cap at the growing plus end, whereas the remaining 
MT consists of GDP-tubulin. The size of the GTP cap 
fluctuates in time and depends on the interplay of GTP- 
tubulin on-rate and hydrolysis rates within the MT. The 
GTP cap stabilizes the MT structure mechanically, and 
catastrophes are triggered by the loss of the GTP cap. 
The catastrophe rate therefore is related to the hydroly- 
sis dynamics within the MT and given by the first pas- 
sage rate to a state with vanishing GTP cap. There are 
different models describing this process [4|-|(|, which will 
be discussed in Sec. Ill Bl in more detail. The dynamic 
instability also limits the ability of MTs to generate poly- 
merization forces Q. 

In a living cell, MTs often cooperate in order to gener- 
ate higher forces. Cooperative MT polymerization forces 
play an important role during mitosis in generating forces 
necessary for chromosome separation [8j. A strong co- 
opcrativity in the dynamics of MTs is also relevant in 
processes regulating the cell length, as has been reported 
in Ref. @ for animal cells. MTs also cooperate in the 
formation of cell protrusions, for example, in neuronal 
growth [1^1 . Therefore, it is important to develop mod- 
els and a theoretical framework for the dynamics of MT 
ensembles under force, which can describe cooperative 
effects in the polymerization dynamics and quantify the 
cooperatively generated MT forces. 

Only recently, has it become possible to study force 



2 



generation by MT ensembles in vitro using MT bundles 
growing against an elastic force, which was realized by 
an optical trap [HI]. Experimental conditions were such 
that no rescue events occurred. The experiments showed 
collective catastrophes of the whole MT bundle: Before 
a collective catastrophe, a large fraction of the MT en- 
semble is growing cooperatively against the elastic force; 
growth often terminates in a catastrophe of the whole 
bundle where all pushing MTs nearly simultaneously un- 
dergo a catastrophe. In the experiments in Ref. [Tlj . 
the maximal polymerization force that is reached before 
a collective catastrophe was measured to grow linearly 
with the number N of MTs in the ensemble. 

In Ref. the experimentally observed collective 

catastrophes, and the linear N dependence, could be re- 
produced in simulations, where polymerizing MTs grow 
against an elastic force and interact by force sharing. A 
theory explaining the characteristic linear dependence of 
the maximal polymerization force on the number iV of 
MTs is, however, lacking. The simulations in Ref. [ll| 
included renucleation of MTs after complete depolymer- 
ization, which led to oscillations in the growth dynamics 
of the MT ensemble. The simulations did not include 
rescue events, which can be relevant in vivo. 

In the present article, we will use a model for the col- 
lective dynamics of N MTs growing against an elastic 
barrier which is very similar to the model that has been 
used in the simulations in Ref. As in Ref. [ll[, the 
ensemble of MTs grow against an elastic force, which 
models the optical traps used in the experiments or the 
elastic cell cortex in vivo. As in Ref. [Hj] , the most impor- 
tant features of the model are that MTs only interact by 
force sharing between the cooperatively pushing leading 
MTs and that the single-MT catastrophe rate increases 
exponentially with its load force. In contrast to Ref. fill ] . 
we will not consider renucleation of shrinking MTs but 
include rescue events into the model, which are an essen- 
tial part of MT dynamics, and study their influence on 
the cooperative MT dynamics. 

Furthermore, we also develop a dynamical mean-field 
theory, which provides a theoretical framework to de- 
scribe the cooperative MT dynamics both in the absence 
and presence of rescue events. It allows us to extract the 
relevant control parameters, such as tubulin on-rate, res- 
cue rate, and MT number, and to investigate their influ- 
ence on dynamics and force generation. It is shown that, 
apart from collective catastrophes, also collective rescue 
events emerge if rescue events are taken into account in 
the single-MT dynamics. The resulting interplay between 
collective catastrophes and collective rescue events gives 
rise to an oscillatory growth dynamics of the entire MT 
ensemble. The dynamical mean-field theory successfully 
describes this oscillatory dynamics as a robust limit cy- 
cle. 

The dynamical mean field theory allows us to calcu- 
late the cooperative polymerization force generated by 
the MT ensemble and its dependence on the MT number 
N, the stiffness of the barrier, and the MT growth pa- 



rameters. Furthermore, we use the dynamical mean field 
theory to investigate the robustness of our results against 
variations of the catastrophe model. This is an important 
question as at least two different catastrophe models have 
been put forward in the literature and have been shown 
to describe experimentally available data on single-MT 
catastrophe rates. Robustness of our results show that 
details of the catastrophe models are not essential for 
force generation by ensembles of MTs but only rather 
general features, such as the exponential increase of the 
catastrophe rate with force, are relevant. We corrobo- 
rate all our mean-field results by microscopic stochastic 
simulations of the full MT ensemble dynamics. 

The paper is structured as follows. In Sec. [Til wc intro- 
duce the model describing the stochastic growth dynam- 
ics of an ensemble of N MTs growing against an elastic 
barrier. We, first, outline the stochastic growth model 
for a single MT, the catastrophe models used through- 
out the work, and, finally, the coupling between the N 
MTs in the ensemble dynamics via force sharing between 
leading MTs. In Sec. IHII we describe our choice of model 
parameters and describe the simulation. 

We then develop the dynamical mean-field theory for 
this stochastic model. In Sec. IIV1 we start with the 
case of zero rescue rate, which has direct applications to 
the experiments by Laan et al. [Til ] and for which the 
cooperative MT dynamics is conceptually simpler to un- 
derstand because there arc only collective catastrophes. 
For zero rescue rate, we discuss the maximal polymeriza- 
tion force and find a linear N dependence for small iV 
or stiff barriers with a crossover to a logarithmic depen- 
dence for large N or soft barriers. Th e results agree with 
the experimental findings of Ref. 
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In Sec. |VJ we introduce the dynamical mean-field the- 
ory for the full problem in the presence of rescue events, 
where the cooperative MT dynamics exhibits both col- 
lective catastrophes and collective rescue events. For 
nonzero rescue rate and a soft barrier, we find stable 
collective catastrophe and rescue oscillations with max- 
imal and average polymerization forces growing linearly 
with N. We show that our theory is applicable to differ- 
ent catastrophe models for single MTs and that our main 
findings are robust for catastrophe models with catastro- 
phe rates exponentially increasing with force. 

Throughout the article, we show that all theoretical 
results are in agreement with microscopic stochastic sim- 
ulations. 



II. MODEL FOR COOPERATIVE DYNAMICS 
OF MT ENSEMBLES 

A. Single MT model 

The dynamic instability causes the MT plus end to 
switch stochastically between a growing (+) and a shrink- 
ing (-) state [13 • In the growing state, GTP-tubulin 
dimers (called monomers in the following) attach and 



3 



detach with rates w on and w ff, respectively, to one of 
the 13 protofilaments. The MT growth velocity in the 
growing state and in the absence of external forces is 
v + = d(ui on — Lu s) [v+ = 1...5 x 10~ 8 m/s], where 
d ~ 8nm/13 is the effective monomer size. 

Under force, the MT growth velocity becomes force 
dependent 



v+(F) 



,-F/Fo 



Wofi 



(1) 



with a characteristic force F = fc^T/d ~ 7 pN. For 
simplicity, we assume that force only affects the on-rate 
of tubulin monomers. Experimental measurements of 
the force-velocity relation in Ref. @] gave a significantly 
smaller value Fq ~ 2 pN. 



B. Dynamic instability and catastrophe models 

The dynamic instability of MTs is associated with the 
loss of the stabilizing GTP cap because of hydrolysis 
within the MT 0. At the catastrophe rate uj c the MT 
loses its GTP-cap in the growing state and switches into 
a state of fast shrinking with a large shrinking velocity 
V- (~ 3 x 10~ 7 m/s). With the rescue rate w r the MT 
switches back from the shrinking into the growing state. 

The catastrophe rate is related to the hydrolysis dy- 
namics within the MT and given by the first passage 
rate to a state with vanishing GTP-cap. The resulting 
catastrophe rate has been discussed based on a model 
for cooperative hydrolysis of GTP-tubulin by Flyvbjcrg 
et al. jj, Q . Similar cooperative models have been pro- 
posed for the hydrolysis dynamics in filamentous actin 
[13l HH . In a cooperative model, hydrolysis proceeds by 
a combination of both random and vectorial mechanisms. 
In a random mechanism, the hydrolysis rate r per length 
(~ 3.7 x H^m-V -1 ) of GTP-tubulin is independent of 
the position of the GTP-monomer in the MT. In a vec- 
torial mechanism, only GTP-monomers are hydrolyzed 
which have already hydrolyzed GDP-monomers at one 
neighboring site. This results in a directed motion of 
GTP- and GDP-tubulin interfaces with mean velocity Vh 
(~ 4.2 x 10~ 9 m/s). The inverse catastrophe rate u)~ x 
can be obtained as mean first passage time to a state 
with a vanishing GTP cap. For a cooperative model, the 
exact analytical result for u> c has been obtained in Refs. 
[3, [1| as implicit function of the growth velocity v + and 
the hydrolysis parameters Vh and r. The exact dimen- 
sionless catastrophe rate a = uj c D^ x ^r^ 2 l z is given by 
the smallest solution of 



Ai'(j 2 - a) = -jAi(-f 2 - a) 



(2) 



with 7 = vD 2 / 3 r 1 ' 3 /2, where v = 14. — Vh and 
D = (v + + v h )d/2 [Ai'(x) = dAi{x)/dx\. Here Ai de- 
notes the first Airy function [l5j]. We use a numerical 
implementation of this exact analytical result for uj c in 
simulations and mean-field calculations: We solve @ to 



calculate the function a = uicD^ 1 ^^ 2 / 3 as a function 
of 7 numerically. From this numerical solution we ob- 
tain the function uj c = uj c (v + ), for which we generate 
an accurate interpolating polynomial of high order. This 
polynomial is used in the simulation to calculate catas- 
trophe rates as a function of (force-dependent) growth 
velocity, v+(F ). The remaining hydrolysis parameters 
Vh and r are fixed during the simulation. 

The catastrophe rate oj c becomes force dependent via 
the force-dependence |T]) of the growth velocity. The re- 
sulting catastrophe rate is a nonlinear and monotonically 
increasing function of the force F (see Fig. [I]), which in- 
creases exponentially above the characteristic force Fq. 
We use the theoretical value F = ksT/d ~ 7 pN for the 
model by Flyvbjerg et al. in the following. 

We will investigate whether collective catastrophe and 
rescue oscillations are robust with respect to the catastro- 
phe model. Different catastrophe models have been pro- 
posed and have been shown to describe experimentally 
available data on single-MT catastrophe rates. One al- 
ternative phenomcnological catastrophe model has been 
proposed by Janson et al. based on experimental data for 
the inverse catastrophe rate, i.e., the average catastrophe 
time r c = 1 /u> c @ . The experimental data show that t c 
increases linearly with the growth velocity v+ such that 
the catastrophe rate is given by 
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bv + 



(3) 



with a ~ 20 s and b : 

Also within the catastrophe model by Janson et a/., 
the catastrophe rate uj c becomes force dependent via the 
force dependence (JTJ) of the growth velocity, and the re- 
sulting catastrophe rate is a nonlinear and increasing 
function of the force F (see Fig. Q]) , which increases ex- 
ponentially above the characteristic force F Q . 

In the simulations of MT ensembles under force in Ref. 
pd| . this catastrophe model is used in combination with 
a characteristic force F ~ 0.8 pN based on the experi- 
mental data in Ref. 0] ; this value is significantly smaller 
than the theoretical value Fq = ksT/d ~ 7 pN. For sim- 
ulations, a 10-fold increased catastrophe rate is used in 
Ref. [llj]. The other parameter values used in Ref. [ll| 
are v on = uj on d = 4.37 x 10~ 8 m/s, which corresponds to 
lj oix ps 70 s -1 @, [ll|. We will use the same set of param- 
eters, apart from the 10-fold increase of the catastrophe 
rate, i.e., a small value Fq ~ 0.8 pN for our simulations 
with the catastrophe model by Janson et al. 

Both the Flyvbjcrg and Janson catastrophe models de- 
scribe available experimental data on single MTs, as has 
been shown in Refs. [3,[f| and [||, respectively. Therefore, 
it is important to investigate whether they also give com- 
patible results for the collective dynamics of MTs under 
force. 

In both models the catastrophe rate uj c decreases as a 
power law over a wide range of growth velocities, 
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(Flyvbjerg), 



UJ C OC V, 



(Janson) (4) 
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FIG. 1: The catastrophe rate uj c (in s" 1 ) as a function of the 
load force F (in units of Fa) for the Flyvbjerg model (with 
Fq ~ 7 pN and for w on = 70 s _1 ), the linear catastrophe 
model (with a = 0.005 s _1 ,& = 8 x 10 5 m _1 , F ~ 7 pN, and 
Won = 70 s _1 ) and the Janson model (with Fo ~ 0.8 pN). 
Inset: The catastrophe rate w c (in s _1 ) as a function of the 
growth velocity i>+ (in ^m/min) for the three different catas- 
trophe models. 



Because of w+ ~ duj on e~ F ^ F ° for large velocities, uj c in- 
creases exponentially with force F above the character- 
istic force Fo in both catastrophe models. This can also 
be seen in the comparison in Fig. [TJ 

Our theory will predict collective catastrophe and res- 
cue oscillations of the MT ensemble for all catastrophe 
models which fulfill two conditions: 

(i) Force dependence via growth velocity: The catas- 
trophe rate uj c = u c (v + ) is a function of the growth 
velocity and becomes force-dependent via the force 
dependence of the growth velocity, uj c = lj c (v+(F)); 
if single-MT catastrophes are related to hydrolysis 
within the MT and the loss of the stabilizing GTP 
cap, oj c should be a decreasing function of v+ and, 
thus, an increasing function of the force F. 

(ii) Exponential force dependence above Fq: The 
resulting force dependence of uj c is such that 
FduJc/dF >• uj c (F) for F > F above the charac- 
teristic force F [see also Eq. This gives rise 
to a catastrophe rate that increases exponentially 
with force above the characteristic force Fo. 

Requirement (ii) is fulfilled for all catastrophe rates de- 
creasing as a power-law uj c oc v^ e (e > 0) with growth 
velocity as in both the Flyvbjerg and Janson models; see 
Eq. ([3]). Accordingly, we will find collective catastrophe 
and rescue oscillations for both models. 

It is possible to consider other types of catastrophe 
models where the catastrophe rate w c = uj c (v+) is a de- 
creasing function of the growth velocity according to re- 
quirement (i) but where requirement (ii) of an exponen- 
tially increase of the catastrophe rate with force is vio- 



lated. One particularly simple example is a catastrophe 
rate which decreases linearly with velocity, 

uj c (v + ) = a - bv + . (5) 

We will demonstrate that collective catastrophe and 
rescue oscillations are indeed absent for such a "linear 
model." 

We will not discuss more elaborate multistep catastro- 
phe models with more than two MT states, which have 
been proposed only recently (Td| . 



C. Model for MT ensemble 

We consider an ensemble of N parallel MTs, directed 
along the x-dircction. The ensemble is growing in a posi- 
tive a;-dircction and pushing against an elastic barrier, as 
shown in Fig. [2] The cooperative dynamics is governed 
by the number n + < TV of leading MTs which push si- 
multaneously in the growing state. 

The elastic barrier is modelled as a spring with equi- 
librium position xq = 1 /jm and a spring constant k in 
the range 10~ 7 N/m (soft) to 10~ 5 N/m (stiff as in the 
optical trap experiments in Ref. [TT| ) - Barrier displace- 
ment by the leading MTs with their tips positioned at 
x > xq causes a force F = F{x) = k(x — xq) resisting 
further growth; for x < xq there is a force-free region. 
We assume that the force F is equally shared between 
all n + leading MTs such that each leading MT is subject 
to a force F/n+. Force-sharing is the only coupling be- 
tween the MTs. In the presence of rescue events, i.e., for 
nonzero rescue rate, we force MTs shrinking to x = to 
undergo rescue. 

Under a shared force F/n+, the growth velocity of 
a MT reduces to v + (F/n + ) = d [w on e- F /" + F ° - w oS ] 
with the characteristic force Fq = ksT/d governing 
monomer attachment; see Eq. ([T]). The force-dependent 
growth velocity also gives rise to a catastrophe rate 
lu c = uj c (v+(F/n + )) increasing with force. All nonlead- 
ing MTs grow with the higher zero force velocity v+(0) in 
their growing state. Therefore, nonlcading MTs, which 
grow force free and fast, "catch up" leading MTs, which 
grow under force with reduced velocity. This mechanism 
supports a state of collective growth, where a relatively 
large number n + of MTs are pushing cooperatively. We 
assume that the shrinking velocity V- is independent of 
force. Because the catastrophe rate depends on force only 
via , the relevant force scale of the problem is set by 
the characteristic force Fo; see (H|). 

This model for the dynamics of the MT ensemble is 
very similar to the model underlying the simulations in 
Ref. [Tll| - In particular, we use the same rules for the 
coupling between MTs by the load force. The most im- 
portant difference is that we include rescue events in the 
single MT dynamics, which have not been considered in 
Ref. [IH. 
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FIG. 2: (Color online) Schematic representation of N — 3 
MTs for n+ = 1. From top to bottom: MTs under force 
F — k(x — xo) grow with velocity v+(F/n+) and unloaded 
MTs with v+(0). MTs shrink with v~ after a catastrophe. 



III. MODEL PARAMETERS AND SIMULATION 

In the simulation model the following parameters for 
MT growth have to be specified: 

(i) The effective GTP-tubulin monomer size d, which 
is d = 8 nm/13 = 0.6 nm for 13 protofilaments. 

(ii) The characteristic force Fq governing the exponen- 
tial decrease (QJ of the single MT growth velocity with 
force. We use F) = ksT/d ~ 7 pN in the Flyvbjerg 
catastrophe model and the measured value Fo ~ 0.8 pN 
in the Janson catastrophe model. 

(iii) We use a GTP-tubulin monomer off-rate of w g = 

6 s _1 [l7| . Then the growth velocity in the absence of 
force v+(0) = d(uj on — w ff) determines the GTP-tubulin 
monomer on-rate cj on (which is proportional to the GTP- 
tubulin concentration GV). 

(iv) Within the cooperative hydrolysis model by Fly- 
vbjerg the catastrophe rate w c is determined by the 
growth velocity v + [according to Q]) with F = ksTjd ~ 

7 pN] and the hydrolysis parameters Vh = 4.2 x 10 -9 m/s 
and r = 3.7 x 10 6 m~ 1 s~ 1 [sj. The exact catastrophe rate 
is calculated from the numerical solution of Eq. ((2|) for 
given u+, Vh, and r, as explained in Sec. Ill B I above. 

(v) Within the hydrolysis model by Janson Q, the 
catastrophe rate uj c is determined by the growth velocity 
v+ [according to (QJ with Fo ~ 0.8 pN] and the parame- 
ters a ~ 20 s and b ~ 1.4 x 10 10 s 2 m- 1 . 

(vi) The shrinking velocity »_=3x 10~ 7 m/s. 
(vh) The rescue rate u) r . 

(viii) For the elastic force F = k(x — xq) on the lead- 
ing MTs we use a spring stiffness k = 10 ~ 7 N/m for a 
soft elastic barrier and k = 10~ 5 N/m for a stiff elastic 
barrier as in the optical trap experiments in [ill ]. The 
rest position xq of the spring is taken as xq = 1 fixa. 
For nonzero rescue rate, we also use reflecting boundary 
conditions at x — 0, i.e., a filament undergoes immediate 
rescue if it shrinks to x — 0. 

The parameters w_, w ff, v^, and r arc fixed in simu- 



lations. We vary the on-rate uj on and, thus, v+(0) and 
the rescue rate u> r within parameter ranges, which are 
selected according to literature values collected in Table 
III 

In the simulation, we integrate the deterministic equa- 
tion of motion for an ensemble of N MTs with continuous 
lengths X{ (i = 1, . . . , N) and include stochastic switch- 
ing between growth and shrinking for each MT. In the 
integration we use a fixed time step At = 0.1 s, which is 
small enough to ensure w CjT .At <C 1. In each time step, we 
have to determine the number n + of leading force-sharing 
MTs. This is done by regarding all growing MTs within 
a distance v + (F/n + )At of the leading MTs as leading 
for the next time step. We can perform two kinds of av- 
erages: Averages (. . .) are taken over many realizations, 
and averages — are timeaverages. 



IV. COLLECTIVE CATASTROPHES AT ZERO 
RESCUE RATE 

We start the analysis with the case of zero rescue rate 
because this case is conceptually simpler to understand 
as rescue events are absent, and there are only collec- 
tive catastrophes to be discussed. Furthermore, this case 
is particularly important because experimental data are 
available: In recent experiments, Laan et al. 11( showed 
that MT ensembles exhibit phases of collective growth 
followed by collective catastrophes, where all leading MT 
nearly simultaneously undergo a catastrophe. The ex- 
periments were performed on short time scales such that 
no rescue events occur. It was also observed that the 
maximal polymerization force before catastrophes grows 
linearly in N. We quantify these features based on a 
dynamical mean-field theory. 

In an ensemble of N MTs the dynamic instability of in- 
dividual MTs leads to stochastic fluctuations in the num- 
ber n + of leading MTs. The force F changes by filament 
growth according to F = kx with x = w+(F/n + ) if the 
ensemble grows (n + > 1) and x = — V- if all MTs shrink 
(n+ =0). In a state of collective growth, a stable mean 
number of MTs arc pushing cooperatively, while the force 
F is increasing by growth against the elastic barrier. If 
the number n + of pushing MTs is reduced by an indi- 
vidual catastrophe, the force on the remaining n + — 1 
leading MTs increases and, thus, their catastrophe rate 
w c (F/n+) increases. A cascade of individual catastro- 
phes — a collective catastrophe — can be initiated until a 
state n + = is reached with all MTs shrinking. This is 
the final absorbing state of the system in the absence of 
rescue events. 

The stochastic dynamics of n + in a growing phase in 
the absence of rescue events is described by a one-step 
master equation with backward rates r n+ = n + ui c (F/n+) 
for decreasing n + by one, which derive from the catas- 
trophe rate of individual MTs under force sharing. In a 
mean-field approach, we replace the stochastic variables 
F and n + by their (time-dependent) mean values (F) and 
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Ref. 


«+(0) [m/s] 


Won [1/s] 


v- [m/s] 


U) r [1/s] 


Drechsel [18] 


(0.7 ... 2) x 10" 8 


(11 ... 32) 


~ 1.8 x 10 -7 


- 


Gildersleeve [19] 


~ 4.2 x 10" 8 


~ 68 


~ 4.2 x 10~ 7 


- 


Walker [20] 


(4 ... 8) x 10 -8 


(63 ... 130) 


~5x 10" 7 


(0.05 ... 0.08)(TUB) 


Laan [11] 


~ 4.2 x 10" 8 


68.25 


- 


- 


Janson [171 


(3 ... 4.3) x 10" 8 


(53 ... 74) 


- 


- 


Pryer [21] 








... 0.5 (TUB) ..0.15 (MAPS) 


Dhamodharan [221 








... 0.07 (Cell) ... 0.085 (MAPS) 


Nakao [23] 








... 0.1 (TUB) 


Shelden [24] 








(0.03 ... 0.2) (Cell) 



TABLE I: Literature values for MT growth parameters v+(0), w on , and uj r . TUB: in vitro results for tubulin solutions, 
Cell: in vivo results, MAPS: effect from MT associated proteins. Values for w on are estimated from measured growth velocities 
via Won ~ v+(0)N/d neglecting ui g. 



(averaging over many realizations of the stochastic 
n + -dynamics) and neglect all higher-order correlations, 
e.g.; set (F/n + ) = (F)/{n + ). In the growing phase, we 
then obtain two coupled mean-field equations, 

d(n + )/dt = -(n + )uj c ((F)/(n + )), (6) 
d(F)/dt = kv + ((F)/(n + )). (7) 

In the mean-field approximation we can calculate the 
maximal polymerization force F meiK (averaged over many 
realizations) that is reached during the mean first passage 
time from n + = N to n + = by solving 

d(F) = d(F)/dt = kv + ((F)/(n + )) 

d(n+) d(n + )/dt (n + )u, c ((F)/(n + )Y [) 

with initial conditions (F) = for (n + ) = N in order to 
find (F) = F max at (n+> « 0. 

Above the characteristic force Fo, the ratio 
v+(F)/lj c (F) decays exponentially because v+(F) 
decreases exponentially and ui c {v+(F)) increases ex- 
ponentially. Therefore, we can solve in two steps: 
(i) As long as the shared force is small compared 
to F , (F)/{n + ) <C Fo, we neglect the force and 
find (F) kv + (0)/ui c (0)ln(N/(n + )). (ii) For 

(F) I (n + ) 3> Fq, on the other hand, the catastro- 
phe frequency increases exponentially, and we can 
assume that d(F) / d{n+) « and (F) remains constant. 

The boundary between regimes (i) and (ii) is deter- 
mined by the condition (F)/(n+) = Fq: Regime (i) 
applies for (n+) > n with n = aW(N/a), where 
W(x) is the Lambert W- function, which is the solution 
of x = We w (for W > —1). The parameter 

a = kv+(0)/u e (fl)F o (9) 

is a dimcnsionless measure for the stiffness of the elastic 
barrier. Because (F) remains constant for < no, 

the resulting maximal polymerization force is given by 

F max = n F = F aW(N/a), (10) 




FIG. 3: (Color online) F max /Fo as a function of N for zero 
rescue rate, w on = 70 s" 1 , and k = 10~ 7 N/m (a ~ 2.7, soft 
barrier) and k = 10~ 7 N/m (a ~ 27, stiff barrier). We com- 
pare results for the catastrophe models of Flyvbjerg (blue, 
solid symbols) and Janson (red, open symbols). Data points 
are simulation results; error bars represent the standard de- 
viation of the stochastic quantity Fmax/Fo. Solid lines are 
numerical solution of the mean field dynamics ([6]) and ([7]). 
Dashed lines are analytical estimates according to (|10|) . 



with a logarithmic asymptotics F max r=s Fqu ln(7V /a) for 
large iV > a or a soft barrier and a quasilinear behavior 
-Fmax ~ FqN(1 — N/a), which is independent of a to 
leading order, for small N <C a or a stiff barrier. The 
mean-field result (fTO]) agrees with numerical solutions of 
the mean field dynamics as given by Eqs. ((6|) and (0 and 
full stochastic simulations both for soft and stiff barriers, 
as can be seen in Fig. [3] 

The parameter a can also be interpreted as a measure 
for the relative speed of the initial and (F) dynam- 
ics according to the mean- field equations ([6]) and (|7|), 
which allows us to give simple arguments for the max- 
imal polymerization force F max : For a <C N (the case 
of a soft barrier), the (n+) dynamics is fast compared to 
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the (F) dynamics. Therefore, (n+) decays approximately 
force free in a time t c ~ 1/cj c (0) In AT from (n + ) = N to 
(n+) = 1. During this time, the force reaches a value 
Fmux ~ fcw + (0)i c ~ FoalnAT. For a ^> N (the case of 
a stiff barrier), the (n+)-dynamics is initially slow com- 
pared to the (F)-dynamics and ~ N until the char- 
acteristic force Fq per MT is reached and the catastro- 
phe rate increases exponentially Up to this point, essen- 
tially N MTs share the force such that F increases up 
t° Fmax ~ FqN until catastrophes set in. This takes a 
time t c ~ NF Q /kv+(0) and A(n+) ~ (n+)i c - A^/a < 1 
is indeed small such that the assumption « ./V is 

consistent. 

In the experiments in Ref . [ll[ , the spring stiffness was 
k ~ 10~ 5 N/m, which gives a ~ 27 such that these ex- 
periments were performed in the quasilinear regime of a 
stiff barrier, where we predict -F ma x ~ FqN for all ex- 
perimentally accessible N (see the upper lines in Fig. [3]). 
This linear increase is in agreement with the experimen- 
tal results but the ratio F max /N is only of the order of 
3 pN experimentally, while Fq ~ 7 pN. This hints at 
a lower value for Fq in the force-polymerization velocity 
relation for MTs; experimentally, a value Fq ~ 2 pN has 
been measured in Ref. Q, which is indeed compatible 
with the experimental results of Ref. • 

In Fig. |21 we compare mean field theory and simula- 
tion results for the Flyvbjerg and the Janson catastrophe 
model for both soft and stiff barriers. For both models, 
we find agreement between mean-field theory and simula- 
tions and, moreover, both models give comparable values 
for generated forces. This demonstrates that results for 
the cooperative force generation at zero rescue rate are 
robust with respect to details of the single-MT catastro- 
phe model. The essential feature entering the mean-field 
theory is the exponential increase of the catastrophe fre- 
quency with force above the characteristic force Fq. 



V. COLLECTIVE CATASTROPHES AND 
COLLECTIVE RESCUE FOR NONZERO 
RESCUE RATE 

We now consider force generation in the presence of 
rescue events. Rescue events were not included in the 
simulation model in Ref. [ll|. Also experiments in Ref. 
[lH were performed on short time scales such that no 
rescue events occurred. However, rescue events are an 
essential part of MT dynamics, and their influence on 
force generation and MT dynamics needs to be addressed. 

In the presence of rescue events, the dynamics will not 
change considerably for A~ <C a, i.e., for a stiff barrier 
because this limit corresponds to a slow (n+) dynam- 
ics, which cannot benefit from additional rescue events. 
Moreover, because the (F) dynamics is fast, rescue will 
not happen before the force- free region x < xq is reached, 
where MTs decouple. 

Therefore, we focus on the influence of rescue events 
for JV > a or a soft barrier corresponding to a fast (n + ) 



dynamics. In this regime, the collective dynamics be- 
comes strongly modified. Apart from collective catastro- 
phes also collective rescue events occur: After a collective 
catastrophe the system is in a state n+ = with all MTs 
shrinking. Individual rescue events lead to n + = 1, but 
a single MT bearing the whole force undergoes an im- 
mediate catastrophe again with high probability. There- 
fore, a cascade of rescue events — a collective rescue — is 
necessary to increase n+ back to a number sufficient to 
maintain stable collective growth. 



A. Simulation results 

Alternating collective catastrophes and collective res- 
cue events give rise to oscillations in the polymerization 
force or, equivalently, the position x of the obstacle. Such 
oscillations with alternating collective catastrophes and 
collective rescue events can be clearly seen in the stochas- 
tic simulation trajectories for the positions x of the MTs 
[see Figs. BJa) and Etc)] and the number n + of leading 
MTs [see Figs.QJb) and0Jd)] as a function of time t. The 
simulation trajectories also show that this phenomenon is 
robust with respect to the catastrophe model and can be 
observed both for the Flyvbjerg and the Janson catas- 
trophe models, which are shown in Fig. @] on the left 
and the right sides, respectively, and exhibit qualitatively 
very similar behavior. Similar oscillations have been ob- 
served in the simulations in Ref. fill ] in the presence of 
MT renucleation instead of MT rescue and for a constant 
force. 

In the simulations we measure the polymerization force 
-Ps.N = (F), which is averaged over time and many real- 
izations as a function of MT number N and of the on-rate 
u; on . This time- averaged polymerization force is also the 
stall force of the MT ensemble. The results are shown 
in Fig. [5] for both the Flyvbjerg catastrophe model [up- 
per row: Figs.[5ja) and[5{b)] and the Janson catastrophe 
model [lower row: Figs. EJd) and[5je)]. The main find- 
ing of the simulations is an approximately linear increase 
of the polymerization force F Si n with the number A" of 
MTs [see Figs.E^a) and[5jd)]. This shows that for large 
MT ensembles, rescue events give rise to much higher 
polymerization forces as compared to the logarithmic A^ 
dependence derived in the previous section in the absence 
of rescue events for a soft barrier (N 3> a). Simulations 
also show an approximately linear increase of the poly- 
merization force with the on-rate w on [see Figs. EJb) and 

We also show numerical results for the time-averaged 
pushing fraction F+ = (n+)/N of MTs as a function 
of the on-rate u> on in Figs. [SJc) and Elf). The push- 
ing fraction increases with on-rate, which demonstrates 
an increasing tendency of MTs to push synchronously at 
higher on-rates, where larger forces are generated. 

Simulation results for the Flyvbjerg model (upper row 
in Fig. [S]) and the Janson model (lower row in Fig. [S]) 
show a very similar linear increase for the polymcriza- 
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FIG. 4: Typical simulation trajectories for N = 10 MTs, 
to on — 70 s~ , and oj r = 0.05 s ; left [(a) and (b)] for the 
Flyvbjerg catastrophe model and right [(c) and (d)] for the 
Janson catastrophe model, [(a) and (c)] Positions of all MTs 
as a function of time t; the obstacle position x(t) is the posi- 
tion of the leading MT. [(b) and (d)] The number n+ of lead- 
ing MTs as a function of time t. Collective catastrophes and 
collective rescue events can be clearly recognized: In a collec- 
tive catastrophe n+ drops to n+ = 1 and x(t) of the leading 
MTs starts to shrink; after a collective rescue n+ starts to 
increase again to values n+ > 1, and x{t) of the leading MTs 
start to grow. 



which are not pushing: The time r denotes the mean 
time that it takes for a MT to rejoin the group of n + 
pushing MTs after undergoing an individual catastrophe 
followed by rescue and force-free growth at a velocity 
v+(0) that is larger than the velocity v+(F/n + ) of the 
leading MTs under force (see Fig. EJ. After a rescue 
time l/ui r the trailing MT has to "catch- up" a distance 
[v+(F/n + )+v-]/uj r to the leading MTs, which kept grow- 
ing with velocity v + (F/n + ). Given a velocity difference 
v + (Q) — v + (F/n + ) to the leading MTs under force, this 
requires a time 

r » w" 1 [1 + (v+(F/n+) + v-)/(v+(0) - v + (F/n+))) 

(11) 

which is larger than the bare rescue time l/uj r . This 
results in a modified mean- field equation for 

d(n+)/dt = -lj c ((F) /(n+)) (n + ) + (r)- 1 (JV - <n+» 

(12) 

where we have to apply the mean held averaging also to 
(r) in Eq. dHJ) : 



v + ((F)/(n+}) + V- 



v + (0)-v + ((F)/(n + )) 



(13) 



Typically (r) is by a factor of 10 larger than the bare 
rescue time l/cj r , 



tion force F Si n with N and very similar results for the 
time-averaged pushing fraction TZjT of MTs, which is in ac- 
cordance with the qualitatively similar simulation trajec- 
tories shown in Fig. |4] for both catastrophe models. This 
further supports that our results are robust with respect 
to the catastrophe model. The absolute values of typical 
forces in Figs. El[a), EJb), EJd), and EJe) and, similarly, 
between typical MT lengths in Figs. |4|a) and EJc) dif- 
fer, however, between the two catastrophe models. The 
reason is that the basic force scale of the problem is the 
characteristic force Fq , above which the catastrophe rate 
increases exponentially, as will be shown below. We have 
chosen the theoretical value Fq = ksT/d ~ 7 pN for the 
Flyvbjerg model and the much smaller value Fq = 0.8 pN 
according to Ref. [ll| with the Janson model. In units of 
the characteristic force Fq, typical forces are very similar 
[see Figs. [SJa), Mb), EJd), and EJe) right scale]. 



B. Dynamical mean field theory 

We will show that all simulation results and the ro- 
bustness with respect to the catastrophe model can be 
explained based on a dynamical mean-held theory. 

In the presence of rescue events, the mean-held equa- 
tion © for (n + ) becomes modified in the growing phase. 
The one-step master equation for n + in a growing phase 
also contains a forward rate g n+ = (N — 7t, + )t _1 for in- 
creasing n + by one. This forward rate is determined by 
a rescue and "catch- up" process for the (N — n + ) MTs, 



C. Limit cycle oscillations and absence of 
bifurcations 

For the further analysis of the mean-held dynamics it 
is advantageous to introduce new variables, the average 
force per MT / and the average fraction v+ of pushing 
MTs, 

f = {F)/N, v + = (n + )/N (14) 

with (F) I (n+) = f /v+- Using these variables, the mean- 
held equations become 

du + /dt = -P + GJ c (f/v+) + (l-V+)/(T) (15) 

df/dt = kv+(f/p+)/N (16) 

We hrst discuss the nullclines of / and v + , i.e., the 
contours in the f-v+ plane along which df/dt = and 
dv + /dt = is satished, respectively. 

The nullclines of / require v+(f = 0, which leads 
to a straight line, 

/ = ^stall (17) 

in the f-v+ plane, where the slope is given by the single- 
MT stall force F sta ii = F ln(w on /w off ) [see Fig. E^a)]. 
The nullclines of v+ are given by 

= g(f, u+) = -v + u c tf/v + )+ 1 . (1 - u+) (18) 
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FIG. 5: (Color online) Upper row: Simulation results using the Flyvbjerg catastrophe model. Lower row: Simulation results 
using the Janson catastrophe model. Both are for a soft barrier (k = 1CP 7 N/m). [(a) and (d)] Simulation results for the 
average polymerization force _F Sj n as a function of the number N of MTs for different on-rates u) on = 30 s- 1 (■), 50 s" 1 (•), and 
70 s _1 (a) at a fixed rescue rate u r = 0.05 s _1 . Solid lines: Mean field estimate F s ,jv = F c /2, see Eq. (|23[) (neglecting -Fmin). 
Dashed lines: Numerical mean field solution including stochastic effects. Black solid line: iV-fold singie-MT stall force F 3 .m = 
NF Bta ,n for Lj on — 70 s _1 . [(b) and (e)] Polymerization force _F s ,n as a function of on-rate oj oi1 for N = 5 (♦), 10 (•), 15 (a), 20 (■), 
and oj r = 0.05 s _1 . Solid lines: Mean field estimate F Sy N = F c (w on )/2. Dashed lines: Numerical mean-field solution including 
stochastic effects, [(c) and (f)] Time-averaged pushing fraction ~ = (n+)/N of MTs as a function of aj on for oj r = 0.05 s" 1 (•), 
0.1 s" 1 (A), and 0.2 s" 1 (■). Solid lines: Solution of Eq. Ip]) for f = f c /2. Dashed lines: Solution of Eq. (T8J for / = F B , N /N 
with _F Si n from the numerical mean-field solution including stochastic effects. Error bars represent the standard deviation of 
the stochastic quantity n+/N. For reasons of clarity we only show error bars for u r = 0.1 s _1 . All other standard deviations 
are of the same magnitude. 
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6: (Color online) Illustration of time scale r from Eq. 



and are independent of N. The shape of the nullclincs 
depends on the functional dependence of the catastrophe 
rate on the force and, thus, on the catastrophe model. 
We will first focus on the Flyvbjerg catastrophe model, 
for which the nullclines of v+ have a characteristic loop 
shape, as shown in Fig. [71 which exhibits two solution 
branches: a stable upper branch v+ M corresponding to 
a collectively growing state with (n+) = Nv +M pushing 
MTs and an unstable lower branch v+. u . For a soft bar- 
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FIG. 7: Nullclines of the mean field equations for (a) the Fly- 
vbjerg model and (b) the Janson model for cj n = 70 s" 1 , 
uj r = 0.05 s" 1 , and k = 10" 7 N/m (soft barrier). The null- 
clines of u+ = (n+)/N axe solutions of Eq. (|18|) and assume a 
loop shape as a function of / = {F)/N. Gray arrows indicate 
the mean-field flow of (n+). The nullcline for / is a straight 
line v+ = ///stall- The critical force is f c — F c /N ~ 3 pN 
for the Flyvbjerg model (a) and f c — 0.4 pN for the Janson 
model (b). Black arrows indicate the stable mean-field limit 
cycle (see text); dashed line indicates effect of stochastic fluc- 
tuations. 
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rier, the (n + ) dynamics is fast, and the force increases 
slowly during collective growth, while v + = v + u is trac- 
ing the stable upper branch of the nullcline. The force 
per MT can increase up to a critical value f c (with a 
corresponding value v c for v + ), where stable and un- 
stable branch join and where the nullcline has vertical 
slope df /dv + = in the f-v + plane. The critical force 
F c = fcN represents the maximal load force for which 
the MT ensemble can maintain a stable state of collec- 
tive growth: for (F) > F c the number (n + ) of pushing 
MTs has to flow spontaneously to a state (n + ) = 0. 

The critical force f c per MT can be obtained from two 
conditions: (i) the nullcline equation g(f,v + ) = 0, i.e., 
Eq. (fT8|) , and (ii) taking a total derivative with respect to 
v+ and using the condition of a vertical slope df /dv+ = 
we arrive at the second condition -^-g(f,v + ) = 0. Be- 
cause U! C (F) increases exponentially above Fo, see Eq. 
((4]), duj c /dF ~ oJc/Fq is a good approximation. The ef- 
fective rescue time (r) = (r)(f /v+) has a much weaker 
force dependence, which we neglect. These approxima- 
tions give 







f 



v+Fo 



'>(//*+) 



V+F 



1 + 



1 



w c (J/ V+ )(t)(J/v+) 



(19) 



It turns out that (for both the Flyvbjerg and the Janson 
catastrophe models) lo c {t) > 1 holds over the entire range 
of forces. In order to estimate f c , we assume lo c (t) 1. 
This leads to an estimate — sw Fo in Eq. (fT9"]) , which can 
be used in the arguments of uj c and (t). Solving the Eqs. 
(US]) and p^|) for f c and v c we find analytical estimates, 



fc 



Fn 



1 



u e {F ) (r)(F ) 
1 

1+lj c (F )(t)(Fq)' 



(20) 
(21) 



According to Eqs. (TTT)) and (|T5|) . the nullclines for / and 
v+ and, thus, the critical values f c and v c are strictly 
independent of N. Therefore, the critical total force 
F c = N f c has to be strictly linear in the number of MTs. 
The critical force is the maximal polymerization force 
that can be generated during polymerization in the pres- 
ence of rescue events. For a soft barrier {N » a), rescues 
thus lead to a significant increase in the maximal poly- 
merization force with a linear N dependence compared to 
the logarithmic dependence derived above in the absence 
of rescue. Moreover, the estimate ([2"0")) for f c predicts 
an increase of the generated force with the on-rate w n 
because this increases i>+ and, thus, reduces uj c and an 
increase with the rescue rate ui r because this decreases 
(t) [see Eq. (HI])]. 

In order to analyze the system for fixed points, 
we compare the lower branch v+j of the nullcline of 
v + with the nullcline (JT7J) of /. The lower branch 
is governed by the exponential increase uj c (F) ~ 



w c (0)exp(cF/F ) with force (with c = 2/3 in the Fly- 
vbjerg and c = 1 in the Janson catastrophe model) 
resulting in w c (0) exp(cf /v+Fq) ~ 1/(t)v + or f/v+ ~ 
Fq ln(l/a; c (0)(T)i/+)/c. This is always at lower forces 
than the nullcline [Eq. [T7[) ] of / because F sta ii = 
F ln(uj on /uj o ff) > F ln(l/w c (0)(r)i/ + )/c. This inequal- 
ity can be violated only at very high rescue rates u> r giv- 
ing rise to a small (r). We obtained that uj r 3> 1/s is 
necessary to obtain a fixed point. Only if this fixed point 
exists and is stable, it can undergo a Hopf bifurcation on 
lowering the rescue rate. We conclude that, for realistic 
parameter values uj t ~ 0.05 s^ 1 , we are always far from 
a Hopf bifurcation. 

The system rather oscillates in a stable limit cycle: Af- 
ter rescue [A in Fig. [7(a)], the pushing force / increases 
with the MT growth velocity because of / = kv + /N, 
while v + = v + , u is tracing the stable branch of the null- 
cline. At the critical force level f c , a collective catastro- 
phe occurs [B in Fig.[7Ja)], where the ensemble is quickly 
driven to collective shrinking with (n + ) = or 1 and 
/ = -kv-/N. 

During shrinking the force level is reduced until an 
individual rescue event can initiate collective rescue at a 
force F m i n [C in Fig. UJa)]. During rescue (n + ) increases 
quickly back to its stable fixed point value [A in Fig.[JJa)] 
closing the limit cycle. 

The collective rescue force F m i n can be calculated from 
the condition that the lower unstable branch of fixed 
points given by Eq. (fT5f intersects the line (n + ) = 1, 
leading to the condition 



N = uj c (F min )(T)(F min ) + l. 



(22) 



Collective rescue typically happens at rather small force 
-Fmin "C Fq such that we find an essentially linear N 
dependence F min ~ N + 0(1). 

The collective mean-field dynamics thus oscillates be- 
tween forces -Fmin and F c . The resulting time- averaged 
polymerization force 



F s ,n = (F) w (F a 



F c )/2 



(23) 



is also linear in N . This is in agreement with the simula- 
tion results [see Figs. [5(a) and [5(c)]. Because F c 3> F mul 
the result F c w NF /(u} c {F ){t)(F )) from Eq. de- 
termines the dependence of the polymerization force -F^n 
on the on-rate uj on and the rescue rate u r . The estimate 
for F c predicts an increase of the generated force with 
the on-rate w on because this increases u + and, thus, re- 
duces uj c . For the velocity dependence ((4]) and assuming 

v + cc uj on (for co on > Won), Eq. ([20)) gives F c oc uH? 
for the Flyvbjerg catastrophe model (and F c oc Lu on for 
the Janson catastrophe model), which is in qualitative 
agreement with the simulation result of an approximately 
linear increase of the polymerization force with the on- 
rate ui on in Figs. [5(a) and [5(d). From the result ([2"0"T) . 
we also predict an increase of the polymerization force 
with the rescue rate uj r because this decreases (r). The 
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pronounced increase of F s ^ with the on-rate w on demon- 
strates that for an MT ensemble, the polymerization 
force can be sensitively regulated by changing the concen- 
tration of available monomers. We also find the collective 
stall force -F Sj n always remains much smaller than the N- 
fold single-MT stall force, F SjN < NF sta u [see Figs. Ufa) 
and EJc)] in contrast to force-sharing filaments without 
dynamic instability, where F Bl n — NF stil ii holds exactly 
|25[. A further confirmation of the mean- field theory 
is provided by simulation results for the time-averaged 
pushing fraction 77^" = (n + )/N in Figs. [S^c) and [5jf ) . 
Mean field results for v + evaluated using the nullcline 
equation (fl8| for / = F s ^/N show good agreement with 
the simulations results. 

The oscillatory limit cycle dynamics, which gives rise 
to collective catastrophe and rescue oscillations, is robust 
against perturbations because the system is far from a 
bifurcation for realistic rescue rates. Only for very high 
rescue rates u r 3> 1/s, does a stable fixed point exist, 
which becomes unstable in a Hopf bifurcation on lowering 
the rescue rate. 

Similar collective catastrophes and rescues are also 
observed in in vitro bulk polymerization experiments 
[26l [27j . In these experiments many MTs synchronously 
polymerize in a solution with GTP-tubulin concentration 
cgtp- All MTs share the available concentration cgtp 
and grow with a velocity v + (cgtp), which decreases if 
GTP-tubulin is consumed. Here, collective catastrophes 
and rescues are caused by sharing the concentration cgtp 
of available GTP-tubulin, resulting in similar collective 
oscillations as force-sharing induces in the present sys- 
tem. 

Finally, we want to note that the collective dynamics 
for N ^> 1 that we described here differ markedly from 
the dynamics of a single MT (N = 1) p}- For a single 
MT rescue does not happen at a particular force level 
-Fmin but after an average time 1 /w r set by the individual 
rescue rate. The resulting N = 1 mean-field equation 
for the average force (F) is v-/tu r = v + ((F)) /u c ({F)) 
Q and equals shrinking and growing distance between 
individual rescue and catastrophe events. 



D. Robustness with respect to catastrophe models 

An essential requirement for the existence of an os- 
cillatory limit cycle is the loop shape of the nullclines 
of according to the stationary mean- field equation 
(fT51) (sec Fig. [7]). Results presented so far have been de- 
rived from the Flyvbjerg model. We obtain a very similar 
loop-shaped nullcline also with the catastrophe model by 
Janson et al. [see Eq. ([3] . The condition for a loop-shape 
nullcline is the existence of a critical force F c , where the 
two solution branches of Eq. (TIgj) merge in a point with 
a vertical tangent. From Eq. (fT5|). we can derive the 
necessary condition 
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FIG. 8: (a) Nullclines of the mean field equations for the 
linear model © with a = 0.005 s -1 and b = 8 x 10 s m" 1 . 
We use u oa = 70 s _1 , uj t = 0.05 s" 1 , and k = 10~ 7 N/m 
(soft barrier). The nullcline for 1/4. = (n+)/JV is not loop- 
shaped. A stable fixed point exists at the intersection with 
the nullcline for /, which is the straight line u+ — //install- 
(b) Simulation trajectory of the position x(t) of all MTs as a 
function of time t for N = 10. 



for the existence of a critical force F c . Therefore, we 
expect the same type of oscillatory limit cycle for collec- 
tive catastrophe and rescue oscillations for a large class of 
catastrophe models which meet the two conditions stated 
in Sec. Ill Bl (i) The catastrophe rate w c = oj c (v + ) is a 
function of the growth velocity only and (ii) the resulting 
force dependence fulfills condition (|24|) . which gives rise 
to a catastrophe rate increasing exponentially with force 
above the characteristic force Fo. Whereas the Flyvbjerg 
and Janson catastrophe models and, more generally, all 
models with u c cx v+ s (e > 0) fulfill condition (|24|) . it is 
violated for the linear catastrophe model Eq. (J5]). 

This explains that the mean-field result of an oscilla- 
tory limit cycle is robust with respect to variations of 
the catastrophe models: We expect qualitatively similar 
behavior for all catastrophe rates oj c (F), which are ex- 
ponentially increasing with force above a characteristic 
force Fo, for example, in the standard catastrophe mod- 
els by Flyvbjerg et al. d, Q or by Janson et al. 0]. 
This explains the robustness observed in the simulation 
results as shown in Figs. 0] and O 

The condition (|24p is violated for the linear catastro- 
phe model Eq. Accordingly, we do not expect to 
find an oscillatory limit cycle with collective catastrophe 
and rescue events. For this type of catastrophe model 
the nullclines are indeed no longer loop shaped, and the 
mean-field theory rather predicts a stable fixed point [see 
Fig. EJa)] . Simulations confirm that collective catastro- 
phe and rescue oscillations are absent for the linear catas- 
trophe model, and we find a rather stationary position 
x of the obstacle and, thus, a stationary polymerization 
force [see Fig. [8jb)]. 



VI. STOCHASTIC FLUCTUATIONS 

The dynamical mean-field theory explains all simula- 
tion results qualitatively. In order to obtain quantitative 
agreement with stochastic simulations, we have to take 
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into account that the maximal force -Fmax for a collective 
catastrophe is typically smaller than the critical-mean 
field force F c [see Figs. Ufa) and 0(c)] because of ad- 
ditional stochastic fluctuations of n + , which reduce the 
first passage time to a shrinking state n+ = [see Fig. 

w- _ 

Improved mean-field results including this effect [28[ 
agree quantitatively with full stochastic simulations [see 
Fig. Ufa) and[5jc)]. We confirm the linear TV dependence 
of the time- averaged force F s n for small TV < 10, and find 
a slightly stronger than linear increase for larger TV. We 
also reproduce the approximately increase of -F s ,n with 
the on-rate w on [see Figs.[5][b) and[5jd)]. 



VII. CONCLUSION 

In cooperative force generation by an ensemble of TV 
MTs, the interplay between force-sharing and MT dy- 
namic instability gives rise to a complex dynamics, which 
can be described in terms of collective catastrophe and 
rescue events. 

We developed a dynamical mean-field theory [sec Eqs. 
(fT5|) and (TIB)) ] which gives a quantitative description of 
the cooperative MT dynamics in terms of two parame- 
ters, the mean force (F) and the mean number of pushing 
MTs (n+), in both the absence and presence of rescue 
events. Using this mean-field theory we identify the rele- 
vant control parameters, such as tubulin on-rate, rescue 
rate, and MT number, and their influence on force gen- 
eration, and we investigate the robustness against vari- 
ations of the catastrophe model. We validated the dy- 
namical mean-field theory by stochastic simulations of 
the MT ensemble dynamics. 

Our main findings are as follows. In the absence of 
rescue events, the maximal polymerization force before 
collective catastrophes grows linearly with TV for small 
TV or a stiff elastic barrier, in agreement with existing 
experimental data ll|, whereas it crosses over to a log- 
arithmic dependence for larger TV or soft elastic barrier 
[see Eq. (jTUJ) and Fig. [3J . This crossover should be ac- 
cessible in experiments by varying the stiffness of optical 
traps. 

In the presence of rescue events and for a soft elastic 
barrier, the dynamics becomes strongly modified: Collec- 
tive catastrophes and rescues lead to an oscillatory stable 
limit cycle dynamics far from a Hopf bifurcation. These 
oscillations should be observable in vitro in experiments 
such as in Ref. [ll| if the MT lengths are sufficient to ob- 
serve rescue events and if the stiffness of optical traps is 
reduced. Moreover, in vivo the behavior of polarized MT 
ensembles can be explored, as has been shown in Ref. Q, 
and our model predicts synchronized growth and shrink- 
age in oscillations if a polarized MT ensemble is growing 
against an elastic barrier such as the cell cortex. 

In the presence of oscillations, we have quantified the 



maximal polymerization force F c = Nf c in Eq. ([20]) and 
the time-averaged polymerization force F s ^ in Eq. (|23p . 
Both forces are linear in TV [see Figs.[5]Ja) and[5]Jd)], and 
the relevant force scale is the force scale i*o, above which 
the MT growth velocity decreases exponentially and the 
MT catastrophe rate increases exponentially. The linear 
TV dependence of forces in the presence of rescue events is 
remarkable because we find an only logarithmic increase 
with TV in the absence of rescue events for soft barriers 
(see Fig. [3]). Nevertheless, even the maximal polymeriza- 
tion force is significantly smaller than the TV-fold singlc- 
MT stall force. This shows that MTs are not optimized 
with respect to force generation because of their dynamic 
instability, even if they cooperate in an ensemble. On the 
other hand, our analysis also shows that force generation 
in MT ensembles is very sensitive to changes of system 
parameters related to the dynamic instability of MTs; in 
particular, it strongly increases with increasing tubulin 
on-rate (and, thus, decreasing catastrophe rate) or in- 
creasing MT rescue rate (see Fig. [5]). The combination of 
both results suggests that a MT ensemble is not efficient 
to generate high forces but that the dynamic instability 
in connection with the ensemble dynamics allows us to 
efficiently regulate force generation through several sys- 
tem parameters. In the living cell, the on-rate can be 
changed by sequestering tubulin-dimers and catastrophe 
and rescue rates are influenced, for example, by micro- 
tubule associated proteins [H, [2l|, [22| ■ 

Our results also have implications for possible mecha- 
nisms which determine the mean length of the MT cy- 
toskeleton, as it has been studied experimentally in Ref. 
Q . For a fixed stiffness k of the opposing elastic force, 
the average force generated by the MT ensemble corre- 
sponds to an average length of MTs. The linear increase 
of the average force with the number TV of MTs suggests 
that the MT length and, eventually, the cell size should 
exhibit a similar linear increase with the number of MTs 
in a polarized MT cytoskeleton if the stiffness of the cell 
cortex remains unchanged. 

In vivo, regulation mechanisms, which will involve the 
kinetics and spatial variation of concentrations of regulat- 
ing proteins, will be relevant for cooperative MT dynam- 
ics and force generation. The present work provides a 
theory to describe the cooperative dynamics arising from 
force sharing and its dependence on various system pa- 
rameters such as tubulin concentration and rescue rates. 
This is a prerequisite in order to explore spatial and tem- 
poral variations of these parameters in regulation mech- 
anisms in future work. 
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